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We study the properties of the Hooke's law correlation energy (E c ), defined as the correlation en- 
ergy when two electrons interact via a harmonic potential in a D-dimensional space. More precisely, 
we investigate the S ground state properties of two model systems: the Moshinsky atom (in which 
the electrons move in a quadratic potential) and the spherium model (in which they move on the 
surface of a sphere). A comparison with their Coulombic counterparts is made, which highlights 
the main differences of the E c in both the weakly and strongly correlated limits. Moreover, we 
show that the Schrodinger equation of the spherium model is exactly solvable for two values of the 
dimension (D = 1 and 3), and that the exact wave function is based on Mathieu functions. 
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I. INTRODUCTION 

Understanding correlation effects remains the central 
problem in theoretical quantum chemistry and physics, 
and the main goal of most of the new theories and mod- 
els in this research area JTJ [2] - In order to gain some 
insight into electron correlation, two-electron model sys- 
tems have always played a key role by shedding new light 
on the relative motion of electrons. 

Two such models are the Hooke's law atom (or har- 
monium, or hookium) [3H6] and the spherium model j7j- 
rB] , In hookium, the two electrons are bound to the 
nucleus by a harmonic potential, while in spherium, the 
position of the electrons are restricted to remain on the 
surface of a sphere. In both cases, the electrons repel 
Coulombically. For these systems, the exact solution 
of the Schrodinger equation can be obtained for some 
discrete values of the confinement parameter (the force 
constant for hookium and the radius of the sphere for 
spherium) [H |5J [TJ . Consequently, these model systems 
(and others [T5rU9] ) have been extensively used to test 
various approximations within density functional theory 
(DFT) [H1I1DH25]. However, the Hartree-Fock (HF) so- 
lution is not always available in closed-form; in hookium 
for example, although accurate solutions have been found 
[27tl28j. no closed-form expression of the ground state HF 
orbital has been obtained yet. 

According to Lowdin the correlation energy (E c ) 
is defined as the error 

E c = E — Ehf, (1) 

which pertains to the HF approximation |30j . E c is a 
function of the external potential V(r), the dimensional- 
ity of the space where the electrons are moving (D), and 
the interelectronic potential w(u), where u = | ti — T"2 1 is 
the distance between the two electrons. 



We define the Coulombic and the Hooke's law cor- 
relation energies as the correlation energies when the 
two electrons interact via a repulsive Coulomb poten- 
tial \w(u) = u^ 1 ] or an attractive harmonic potential 
[w(u) = oj 2 u 2 /2], respectively. 

We have recently shown that the Coulombic correla- 
tion energy is rather insensitive to V{r) (at least in the 
high-density limit |31j). but strongly dependent on the 
number of degrees of freedom of the electron pair (£>). 
However, the question of how the interelectronic poten- 
tial w(u) influences the correlation energy has not yet 
been addressed, and this is indeed the purpose of this 
article. 

In the following, we will consider the l S ground state 
of the Hooke's law analog of D-hookium (known as the 
Moshinsky atom [32] ), and Z?-sphcrium (labeled as D- 
HL-sphcrium in the following). We adopt the convention 
that a D-sphere is the surface of a (D + l)-dimensional 
ball. Thus, in 2-spherium, for example, this is the sur- 
face of a three-dimensional ball. For D = 3, both the ex- 
act and HF solutions of the Moshinsky atom are known 
|32) . We will generalize these results for any value of 
D in Section [TTJ and derive the two-electron probabil- 
ity distributions in both the position and momentum 
space. In Section |III[ we analyze the energy behavior 
of D-HL-spherium in both the weakly and strongly cor- 
related regime. Moreover, we demonstrate that the exact 
solution of the Schrodinger equation can be found for two 
values of the dimensionality (D = 1 and D = 3). Atomic 
units are used throughout. 

II. MOSHINSKY ATOM 
A. Exact solution 

The Moshinsky atom (32| is defined by the Hamiltonian 
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FIG. 1. — in the Moshinsky atom as a function of lo (solid). 
The small-cj (dashed) and large-o; (dotted) expansion are also 
represented. 




co (a.u.) 



where lo 2 is the force constant between the two electrons. 
The Hamiltonian ([2| can be separated into the extracule 
and intracule coordinates, which reads respectively 



A = 



n +r 2 

V2 ' 



A = 



n - r 2 

V2 ' 



yielding 



ff=i[-Vi + A 2 ]+i[-V| + (2^ + l)A 2 ] 



(3) 



(4) 



For S states in a D-dimensional space, the Laplace oper- 
ator is given by [S3] 



V 2 = 



c) 2 



D — 1 d 



dr 2 r dr 
This leads to the exact wave function 



V>(a,a; 

and energy 



_ (2^ 2 + l)^/^_ u ^_i^ TTA2 



r D/2 



E 



D 



(5) 



(6) 



(7) 



Eq. |7]) reveals that the exact energy is linear with re- 
spect to the dimensionality D 7 due to the separability of 
the Hamiltonian pi). 



B. Hartree-Fock approximation 

According to the Hartree-Fock (HF) approximation 
[35], the total HF wave function of the singlet state is 
defined as 



ipUF{ri, r 2 ) = </?hf Oi) mv (r 2 ) , 



(8) 



where ^hf is the HF orbital, eigenfunction of the Fock 
operator 



F=-Wi 



1 



^{r 2 )rldr 2 , (9) 



associated with the eigenvalue 

D 3lo 2 + 2 



4 V^ 2 + 1 ' 



It can be shown that 



^hf(t-) = ^ D/A 
Then, ^ is easily recast as 



_ ( a ,2 + 1)J /8^, v _ Ir2 



i/W(A, A) 



(^ + 1 ) g/4 r ^7^(A 2 + ^) 



r Z5/2 



(10) 



(11) 



(12) 



The total HF energy, which also behaves linearly with D, 
is 



Eh 



HF 



Wi 



(13) 



C. Correlation energy 



According to ([I]) , the explicit expression of the corre- 
lation energy of the Moshinsky atom reads 
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(l - 2y/u 2 + 1 + ^2lo 2 + 1) , (14) 



which obviously decreases linearly with D. 

Eqs. ([7]), ( fl~3| ) and (14 1 yield the small-w expansion 



= 1 
= 1 



E 
D 

Ehf 
~D~ 

E c _ lo 4 
~D ~ ~~8 4 

and the large-w expansion 
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i?HF 1 



2w 
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- 1 \L) 



n/2-4 
8w 



(15) 
(16) 
(17) 

(18) 
(19) 
(20) 



The correlation energy of the Moshinsky atom, repre- 
sented in Figure [TJ is quartic for small lo (weakly corre- 
lated regime) , and it decreases as 00 in the large-o; regime 
(strongly correlated regime). In comparison, the correla- 
tion energy of _D-hookium in the weakly correlated limit 
tends to a constant [5] [3T] . 



3 



FIG. 2. V (solid), Phf (dashed) and AT 3 (dotted) in the 
Moshinsky atom as a function of u for the unit force constant 
and D = 3 
!P(u) (a.u.) 




FIG. 3. M (solid), TMhf (dashed) and AM (dotted) in the 
Moshinsky atom as a function of v for the unit force constant 
and D — 3. 
M(v) (a.u.) 




In the strongly correlated regime (w — > oo), the exact 
and HF wave functions of the Moshinsky atom become 



V>(A,A) = 
V>hf(A, A) = 



(2w) 



Z5/8 



-e 2 e 75 



associated with the energies 



wjD/2 -i^(A 2 + A 2 ) 



E _ uj 1 

E'hf 1 
^ = " + 2^ 



(21) 
(22) 

(23) 
(24) 



In (23), the second term (1/2) is related with the motion 
of the center of mass, while the first term (w/%/2) is as- 
sociated with the zero-point oscillations of the electrons. 
Indeed, in the strongly correlated regime, the electrons 
oscillate with an angular frequency D^/2uj. Such phe- 
nomena are ubiquitous in strongly correlated systems, as 
demonstrated by Gori-Giorgi, Seidl and their coworkers 

[HI23IM1I3M]. 

As one can see, the HF solution does not describe prop- 
erly these oscillations, and thus exhibits a wrong behavior 
in the large-w limit. 



D. Position Intracule 

To study further the correlation effects in Hooke's law 
systems, we have determined the position intracule V(u) 
of the Moshinsky atom (D > 2). V(u) gives the proba- 
bility of finding two electrons separated by a distance u 
[UJ |32], and is defined by 



where £l u is the angular component of u, and S the Dirac 
delta function [33]. From ^ and (j8j, it follows 



V{u,u>) 
7>hf(w,w) 



2 f-ir(#; 



(26) 



2f-F(f) ■ ^ 



where T is the gamma function 34 . Surprisingly, we 
have the relation 



(28) 



which means the HF intracule is related to the exact 
one by a scaling of the confinement strength. This kind 
of relation could be useful in the future development of 
intracule functional theory [43TI48] . 

In analogy with the Coulomb hole [3T], we define the 
Hooke hole as the difference between the position intrac- 
ule obtained from the exact wavefunction and the corre- 
sponding HF one: 



AP(u,w) =V(u,u)-V kf (u,l>). 



(29) 



Contrary to the Coulomb correlation [JT] 03] , in the 
Moshinsky atom the correlation increases the likelihood 
of finding the two electrons close together and decreases 
the probability of larger values of u. It implies that the 
Hooke hole is positive for small u and negative for larger 
u. This behavior is due to the attractive nature of the 
harmonic potential between the electrons. It is illustrated 
in Figure [2] where we have reported the position intrac- 
ules and the Hooke hole of the Moshinsky atom for to = 1 
and D = 3. 



E. Momentum Intracule 



V(u)= ( ^(n.raj^dri-ral-ttjdnt,, (25) 



The momentum intracule [50] . which gives the proba- 
bility of finding the two electrons moving with a relative 



4 



momenta v, has been computed to get some insight into 
the relative momenta of the electrons: 



M(v) 



(Pi,P2) 2 s{\pi-p 2 \-v)dn v , (30) 



(n„) 



where 4>(pi,P2) is the momentum wave function, pi and 
Pi are the momenta of electrons 1 and 2, and Cl v is the 
angular component of v. 4>(pi,p 2 ) is obtained from the 
position wave function by a Fourier transform: 



0(A,A) 

0hf(A, A) 
where 

A 



(2co 2 



-D/8 



T D/2 



-iA 2 / ■, 

2 Ji g 2V2ui 2 + l 



T D/2 



-l v / 



(31) 
(32) 



(33) 



Pi +P2 z Pi -P2 

are the extraculc and intracule coordinates in momentum 



space. Using (31 1 and (32 1, the exact and HF momentum 



intracules are given by 
M(v,uj) = 



. 2 / 



Mkf{v, uj) 



2T-in 
(co 2 + iy D/4 v D - 1 -- 



(34) 



(35) 



From (34) and (35), we show that the property (|28[) re- 



lated to the position intracule is still valid for the mo- 
mentum intracule: 

Mhf(«,w) =M(v,^=). (36) 
v2 

The results are gathered in Figure [3j where we have rep- 
resented the momentum intracules and the Hooke hole in 
momentum space 

AM(v,u) =M(v,uj) -Mhf(«,w). (37) 

The exact and HF momentum intracules exhibit a simi- 
lar shape to the position intracules, but the Hooke hole 
in momentum space is drastically different from its po- 
sition space counterpart. Indeed, the correlation in mo- 
mentum space favors electrons moving with high relative 
momentum, while correlation in position space evidences 
a decrease of the probability of finding electrons further 
apart. Once again, this behavior is different in hookium, 
in which correlation favors both lower and higher relative 
momenta [27] ■ This is also due to the attractive nature 
of the interelectronic potential. 

III. HOOKE'S LAW SPHERIUM 
A. Expansion of the exact energy 

In terms of the interelectronic angle (0 < 9 < it), the 
Hamiltonian of two electrons on a D-sphere of radius R, 



and interacting with a force constant fl 2 is [T2HT1] 

d 



H = 



1 

R 2 



de 2 



(D- 1) cot6>— + frir (1- cos 6). 
dO 

(38) 

After the energy scaling E 4- R 2 E, and the definition 
of the dimensionless variable w = flR 2 , The Hamiltonian 



( 38 ) reduces to 



H 



de 2 



+ (D- l)cot0 



dO 



+ co 2 (l-cos0). (39) 



Then, perturbation theory [1 m CE3 GIO ED [52] can 

be applied, and expanding E up to the second-order in 
uj 2 , it is straightforward to show that 



E 



D(D + 1) 



(40) 



which is valid for the small-o; regime (weakly correlated 
limit) . 

For large uj (strongly correlated limit), the potential 
dominates the kinetic energy and the electrons tend to lo- 
calize on the same side of the sphere, and oscillate around 
their equilibrium position (zero-point oscillations) . How- 
ever, in Coulombic systems, the localization takes place 
on opposite sides of the sphere (formation of a Wigner 
molecule [S3]), because of the repulsive nature of the 
Coulombic interaction. 



In this limit, the Hamiltonian (39) becomes, for small 
oscillations (0 ~ 0): 



dff 2 



D - 1 d 
6 d9 



where we have used the Taylor expansions: 



n 1 
cot^ — 



cos 9 



9 2 

1 

2 



(41) 



(42) 



The corresponding ground state eigenfunction and eigen- 
value of (41 ) are 



^°\9) 



1 



2-- 1 r(f) 



E (P) 



Duj 



(43) 



The electrons are localized on the same side of the sphere, 
and oscillate around their equilibrium position with an 
angular frequency of D\[2uj. One notes that the electrons 
oscillate with the same angular frequency in both the 
Moshinsky atom and D-HL-spherium (see above). 
The first-order correction of the kinetic energy 



D-l d 
3 d9 



(44) 



and the zeroth-order wave function ( 43 1 yield the asymp- 
totic expansion [12l [13] 



E 



Duj D(D-l) 
7^ 6 



(45) 
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FIG. 4. ip„ as a function of 9 for the unit force constant, 
and various n [n=0 (solid), 1 (dashed), 2 (dotted) and 3 (dot- 
dashed)]. 

(a.u.) 
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(b)D- 



which coincide with the Mathieu differential equation [MJ 
[55] for D = 1 and 3, due to the cancellation of the second 
term in brackets in (48). 

Then, the cases D = 1 and D = 3 are exactly solvable, 
and the exact wave functions of the nth excited state are 



^ D (0) = Se (&a(n-i),-2w 2 ,0 , 



'2(n-l)) 



2w 2 



(49) 
(50) 



where Se is the the odd Mathieu function, and b 



2(n-i) 



its 



characteristic value (n £ N is the number of nodes in the 
wave function between and n) [M] . The energies of the 
nth excited state are given by 



52(n-l) 



? 2(n-l) 



-1, 



(51) 
(52) 



The Taylor expansions of ( 51 ) and ( 52 ) demonstrate that 
the ground state energy behaves consistently with Eqs. 
(|40b and (45) for small- and large-w, respectively. Figure 



141 shows the ground state and the first three excited states 
for D = 1 and D = 3 (u) = 1 in both cases). 



C. Hartree-Fock approximation 



which shows that the next term of this expansion is dif- 
ferent in the Moshinsky atom and Z3-HL-spherium than 
their Coulombic analogues. 

For comparison, in D-spherium (the Coulombic ana- 
log of D-HL-spherium) , one can eventually show that 
the angular frequency is equal to D/(2R 3 / 2 ), and the 
asymptotic expansion reads 



E 



1 



2R 4i? 3 / 2 



D D{9D - 14) 



MR 2 



(46) 



where the first term in ( 46 ) represents the classical me- 



chanical energy of two electrons sitting on opposite sides 
of a sphere of radius R, and the third term is the first 
anharmonic correction. 



B. Exact solvability for D — 1 and D 

Using the ansatz 

ip(0) = (sinfl) 15 ^ x(6») 



(47) 



Eq. ( 39 ) is recast as 
d 2 X (0) 



d9 2 



D - 1 
2 



(D-l)(X»-3) 



cot 2 9 



Following our previous work [T31 [3T] , it is straightfor- 
ward to show that, for D > 2, the HF wave function and 
energy are 

^hf(^) = V D 2 +1 J , E HF = uj 2 , (53) 



D + l ' 

27T 2 



which yields an uniform electron density over the surface 
of the hypersphere. 



D. Correlation energy 

Following the results of the previous section, for small- 
er , the correlation energy for two electrons on a D-sphere 
behaves as 



E r 



D(D + 1) 



(54) 



which, like the Moshinsky atom, is quartic in w, but ex- 
hibits a different behavior with respect to D. The corre- 
lation energy in 3-HL-spherium is represented in Figure 

E 



X (0)-u 2 (1 - cos 



From Eq. (|54|), it is obvious that E c behaves quadrat- 
ic * 
How- 



isalWNwi^h /respect to D in the large- D limit [E c 
l£fl&f, W<i 5-spherium [E c ~ -1/{8D 2 )] [3T]. H. 



(48) ever, the prefactor is different. 
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FIG. 5. — E c in 3-HL-spherium as a function of uj (solid 
curve). The small- and large-o; expansions are also repre- 
sented (dashed and dotted curves, respectively). 

-E c (a.u.) 




IV. CONCLUSION 

In this article, we have studied the Hooke's law corre- 
lation energy of two model systems: the Moshinsky atom 
and the spherium model. 

We have shown that, in the weakly correlated regime 
(small-w limit), the correlation energy is quartic in u in 
both systems, but behaves differently with respect to the 
dimensionality of the space. This feature reveals the dif- 



ference between the Coulombic and the Hooke's law sys- 
tem. Indeed, the correlation energy of both D-hookium 
and D-spherium tends to a constant in the high-density 
limit. 

In the strongly correlated regime (large-cj limit) [3H 
HO] , the leading terms of the asymptotic expansion in the 
Moshinsky atom and D-HL-sphcrium are identical, and 
they represent the zero-point oscillations of the electrons 
when the kinetic energy tends to zero. This could be 
viewed as an "attractive" version of the Wigner crystal- 
lization [53], which involves electrons in the low-density 
limit (quantum dots). 

Moreover, we have shown that the Schrodinger equa- 
tion of the D-HL-spherium reduces to a Mathieu differ- 
ential equation [Ml [55] for two specific values of the di- 
mension: the model of two electrons on the surface of a 
sphere and interacting via a hookean potential is exactly 
solvable for D = 1 and 3, and the exact wave function is 
based on Mathieu functions. 
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